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Membranes  composed  of  am- 
phiphilic  molecules,  such  as  the 
monolayers  of  surfactant  mole¬ 
cules  at  oil-water  interfaces  in  microemul¬ 
sions,  the  lipid  bilayers  that  form  biological 
membranes,  as  well  as  the  layers  of  surfac¬ 
tant  molecules  in  recently  studied  lyotropic 
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liquid  crystals,  are  highly  flexible  (nearly 
tensionless)  surfaces.  Membranes  play  a  cen¬ 
tral  role  in  determining  the  architecture  of 
biological  systems  and  provide  the  basic 
structural  element  for  complex  fluids  such  as 
microemulsions;  an  understanding  of  the 
statistical  mechanics  of  these  .self-avoiding 
surfaces  is  therefore  of  considerable  impor¬ 
tance  (I,  2). 

In  most  cases  of  interest,  these  mem¬ 
branes  are  fluid,  which  means  that  the  mol¬ 
ecules  can  diffuse  rapidly  within  the  mem¬ 
brane  surface  and  possc.ss  no  reference 
lattice.  In  the  absence  of  a  lateral  tension  (J, 


4),  the  shajie  of  the  membrane  is  governed 
by  its  bending  rigidits'  k.  A  membrane  of 
linear  size  L  exhibits  transverse  fluctuations 
(.^)  of  extension  /.  j  ~  {k'l’iK)'  ^l.  (where  k 
is  Boltzmann’s  constant  and  7'  is  tempera¬ 
ture)  on  length  scales  small  compared  to  the 
persi.stence  length  (6)  4,.  =  *'  expliK/l'/'), 
where  a  is  a  short-distance  cutoff  and  i  = 
4tr/3  (7,  S).  Htiwever,  as  L  approaches 
shape  fluctuations  have  been  predicted  (7.  S) 
to  reduce  the  bending  rigidits'  and  lead  to  a 
renormalized  rigidirv'  =  k  -  {kT/c) 
In(Llii)  in  the  limit  of  small  I’V'/k.  On  length 
scales  L  ~  the  membrane  should  have  an 
effective  bending  rigidity  of  the  order  of  feT. 
At  larger  length  scales  (L  >>  mem¬ 
branes  fluctuating  at  constant  area  are  ex¬ 
pected  to  have  an  extremely  small  bending 
rigidity'  and  behave  as  crumpled  (9)  objects 
characterized  by  the  absence  t)f  long-range 
orientational  order  of  normals  erected  per¬ 
pendicular  to  the  l(x;al  surface  elements.  For 
a  fixed  topology,  the  .scaling  behavior  of 
self-avoiding  fluid  membranes  at  these 
length  scales  is  expected  to  be  the  same  as 
that  of  a  branched  ptilymer  (2,  10,  tt). 

Several  aspects  of  this  scenario  have,  how¬ 
ever,  '  ot  yet  been  verified.  Indeed,  there  is 
link  evidence  that  self-avoiding  fluid  mem¬ 
branes  with  a  finite  bare  bending  rigidity 
really  do  crumple.  It  has  been  speculated  (2) 
that  self-avoidance  may  stabilize  the  effective 
bending  rigidity  at  some  finite  value  k  —  kT 
and  therefore  prevent  crumpling.  Further¬ 
more,  there  is  .some  controversv  ctmeerning 
the  universality  classes  of  the  various  models 
used  to  describe  self-avoiding  random  sur¬ 
faces  (with  K  =  0).  The  most  widely  studied 
models  for  these  surfaces  are  constructed  by 
taking  the  elementary  2-cells  (plaqucnes)  on 
a  hypercubic  lanice  and  gluing  them  togeth¬ 
er  in  such  a  way  that  each  edge  is  shared  by 
exactly  two  plaqucnes  (/(?).  For  fixed  topol¬ 
ogy,  the  long-length-scalc  behavior  of  this 
cla.s.s  of  surface  has  bexu  shown  to  be  that  of 
a  branched  polymer:  the  radius  of  gy'ration 
of  a  surface  r)f  area  S  scales  as  ~  S'’ 
(72),  with  V  =  1  in  spatial  dimension  <7  =  3 
(70).  Recently,  however,  randomly  triangu¬ 
lated  surfaces  of  the  type  we  consider  here 
have  been  invc.stigatcd  for  k  =  0  with  the 
use  of  Monte  Carlo  techniques  {13).  A  value 
V  =  0.8  was  reported,  implying  that  such 
surfaces  belong  to  a  different  universality 
class.  Given  the  generality  of  the  cntropic 
mechanism  behind  branched  polymer  be¬ 
havior,  this  discrepancy  is  troubling  (70, 
77). 

In  this  report  we  present  evidence  that 
fluid  membranes  arc  crumpled  at  sufficiently 
large  length  scales  for  any  value  of  the  bare 
bending  rigidity  and  that  the  crumpled  state 
docs  indeed  exhibit  branched  polvmer  be¬ 
havior.  Our  conclusions  arc  based  on  exten- 
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Rg.1  .  Mean-squared  radius  of  gvration  versus 
the  number  of  mtrnomers  N  frrr  \  =  0;  (O)  open 
membrane,  ff,  =  2.0,  frec-edge  boundars'  condi¬ 
tion;  and  ( X )  vesicle,  fj;  =  2.8.  The  .solid  line  is  a 
plot  of  Rg  —  N. 


sive  simulations  of  a  simple  molecular  model 
for  self-avoiding  fluid  membrane.s  for  a  wide 
range  of  bare  bending  rigidities  The  simu¬ 
lations  were  carried  out  with  a  simple  string- 
and-bead  model  for  randomly  triangulated 
two-dimensional  surfaces  embedded  in  three 
dimensions  (13-17).  Both  planar  (with  free- 
edge  boundary  conditions)  and  spherical 
(vesicle)  topologies  were  considered.  In 
both  ca.ses  the  surface  is  modeled  by  a 
triangular  network  of  N  hard-sphere  parti¬ 
cles  of  diameter  tr  =  1 .  The  energy  assigned 
to  a  particular  configuration  is 

P3€=  2)  ("a'n-,  -  1)  + 

(o.7> 

2;  Mr,  -  rfl)  +  21  l^Hcdr,  -  r,|)  (1) 

(I'u) 

where  3  =  l/kT.  In  Eq.  1,  the  first  sum  runs 
over  pairs  {a,y}  of  unit  vectors  {n^J  erected 
perpendicular  to  each  elementary  triangle  in 
the  lattice.  The  second  summation  is  over 
neighboring  pairs  (ij)  of  atoms  (located  at  r, 
and  tj)  in  the  array  interacting  through  a 
tethering  potential  V,(r)  that  vanishes  for  0 
<  r  <  and  is  infinite  otherwise.  Finally, 
the  third  summation  is  over  all  pairs  of 
atoms  interacting  through  a  repulsive  hard 
core  potential  Vm-ir)  that  vanishes  for  r  >  1 
and  is  infinite  otherwise.  Self-avoidance  is 
enforced  by  choosing  an  fj,  <  fo'**;  the 
value  of  depends  on  the  updating 
procedure  that  is  implemented  (see  below). 
The  continuum  limit  of  the  bending  energy 
term  in  Eq.  1  can  be  shown  to  be 

1  f 

Eb  =  -  K  (H^  -  2K)  (2) 

where  dS  denotes  the  surface  element  and  f\, 
is  just  Helfrich’s  curvature  elastic  energy 
with  bending  rigidity  k  =  \/V3  (18,  19) 
and  Gaussian  rigidity  k<.;  =  -k;  H  =  r ,  -f  Cj 
and  K  =  f,f2i  where  r,  and  f2  are  the 
principal  curvatures. 


Our  Monte  Carlo  ujidating  procedure 
consists  of  two  steps:  First,  we  attempt  to 
sequentially  update  ibc  jv-ciriov  vector  of 
each  monomer  bv  a  random  increment  in 
the  cube  |  -s..'  )',  where  .<  is  a  constant  to  be 
chtjsen  below.  The  probabilistic  decision 
whether  to  accept  the  move  is  made  by 
comparing  the  initial  and  final  energies  of 
the  system.  We  chose  s  so  that  ~50%  of  the 
updating  attempts  were  accepted.  Second, 
we  attempt  to  flip  N  randontlv  chosen 
bonds.  A  bond  flip  consists  of  deleting  a 
tether  and  constructing  a  new  one  between 
the  two  previously  unconnected  vertices  of 
the  two  adjacent  triangles.  The  flip  is  accept¬ 
ed  with  the  probability  given  bv  the  Boltz¬ 
mann  factor  if  all  vertices  have  a  minimum 
of  three  neighbors  (20).  Note  that  this  pro¬ 
cedure  docs  not  flip  boundary  tethers,  so 
that  the  open  membranes  we  simulate  have 
perimeters  of  coastant  length.  A  more  de¬ 
tailed  discassion  of  the  bond-flipping  proce¬ 
dure  can  be  found  in  (13  17). 

In  order  to  iasure  reparametization  in¬ 
variance  we  use  a  discretized  version  of  the 
invariant  measure  /D[r]  =  fl^  [/dr(^Je' ■*- 
(4)1  in  our  simulations,  where  r(4)  is  the 
three-dimensional  coordinate  vector  of  a 
point  on  the  surface  with  internal  coordinate 
4  and  g  is  determinant  of  the  metric  tensor 
(21).  We  also  choose  a  metric  in  which  all 
internal  lengths  are  equal,  that  is,  all  trian¬ 
gles  on  the  surface  are  equilateral  and  of  area 
one  so  that  the  volume  o),  of  the  dual  image 
of  vertex  /  (which  is  the  discrete  analog  of 
the  invariant  volume  /  d^^\Q)  is  propor¬ 
tional  to  the  coordination  number  of  this 
vertex:  u),  =  qj$.  On  the  triangulated  surface 
the  curvature  is  concentrated  at  the  vertices; 
with  the  present  normalization  the  scalar 
curvature  Rf  =  ir(6  -  qi)lq,  and  the  discre¬ 
tization  of  the  mea.sure  is  fl,  di^  ( 14-17). 

Simulations  were  performed  w'ith  both 


=  \  2.0  and  \  2.X.  In  the  former  case  we 
chose  s  =  0. 1  and  in  the  latter  .v  =  0.15.  This 
choice  of  p,ir,inietcrs  ensures  self-avoidance 
when  updating  the  monomer  ccxirdinates. 
However,  there  is  only  weak  self-avoidance 
during  Ixind  flipping  for  f„  =  V2.8.  In¬ 
deed,  it  is  easy  to  show  that  implementing 
complete  self-avoidance  would  require  tak¬ 
ing  f„  at  least  less  than  V  8;5  in  this  case.  We 
have  ill  fact  observed  vesicle  inversion  for 
=  V2.8  when  \  =  0.  Nevertheless,  data  for 
the  mean-squared  radius  of  gyration  for 
zero  bending  energy  indicate  that  this  diKs 
not  influence  the  asymptotic  fractal  dimen¬ 
sion  of  the  crumpled  phase.  Furthermore,  it 
docs  not  qualitatively  att'ect  the  behavior  for 
finite  X.  The  primary'  advantage  of  using  a 
larger  tethering  length  is  that  the  dilfusivity 
of  the  monomers  in  the  membrane  surf  ace  is 
substantially  enhanced:  the  bond-flip  ac¬ 
ceptance  rate  is  a  factor  of  2  greater  tor  f  ,, 
=  VT8  than  for  \/2d)  when  X  =  0;  the 
effect  is  greater  for  larger  values  of  X. 

In  Fig.  1  we  plot  data  ftir  R~^  versus  the 
number  of  monomers  for  the  case  of  zero 
bending  rigidity,  X  =  0.  Averages  were 
taken  over  1  x  10^  to  2  x  10^  Monte  Carlo 
steps  per  monomer.  Data  are  shown  both 
for  vesicles  with  f,,  =  planar 

membrane.s  with  f„  =  V-^.O.  For  both  to¬ 
pologies,  the  data  are  consistent  with  R^  ~ 
N''  (12)  with  V  =  1.  Furthermore,  typical 
configurations,  such  as  the  one  illustrated  in 
Fig.  2,  indicate  a  “branched  polymer”  stnic- 
turc  in  which  there  are  long  arms  extending 
out  in  various  directions.  This  result  sug¬ 
gests  an  entropic  mechanism  that  favors 
treelike  ramified  objects,  which  iS  in  agree¬ 
ment  with  results  obtained  for  random  sur¬ 
faces  constructed  of  elementary  plaqucites 
on  a  cubic  lattice  (10).  For  vesicles,  our 
results  are  consistent  w'ith  S  —  R^  and  I’  — 
(w'here  S  and  Fare  the  surface  area  and 
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Fig.  3.  (A)  Specific  heat  obtained  by  using  Hq.  3 
tor  vesicles  with  =  \/J^  and  /V  =  (□)  47,  ( x ) 
127,  and  (  +  )  247  monomers.  (B)  Eft'ective  bend¬ 
ing  rigidity  i<  =  obtained  by  averaging  over 
V,  =  (N  -  l)/2  contiguous  monomers  on  the 
vesicle  surface  by  using  Eq,  5. 

volume  enclosed  by  the  vesicle). 

Fluctuations  are  suppressed  for  finite  \, 
and  as  seen  in  Fig.  2,  there  appears  to  be  a 
cros,sover  from  crumpled  behavior  for  \  <  < 
1  to  a  rough  but  extended  phase  for  \  >  >  1 . 
Further  evidence  for  this  freezing  out  of 
degrees  of  freedom  is  given  by  the  behavior 
of  the  specific  heat  per  monomer 

C  =  -[<£•')-<£)']  (3) 

where  E  =  n„'np.  A  plot  of  data 

obtained  for  vesicles  with  fp  =  con¬ 

taining  N  =  47,  127,  and  247  monomers 
(Fig.  3A)  .shows  that  rhere  is  a  pronounced 
peak  near  k,.  =  1.25.  Such  a  peak  could 
signal  the  existence  of  a  phase  transition 
between  a  crumpled  phase  for  low  bending 
rigidities  and  a  flat  phase  for  large  ones.  In 
fact,  it  is  argued  in  ( 17)  that  such  a  transition 
does  indeed  <x:cur  in  membranes  without 
self-avoidance.  Flowever,  the  height  of  the 
peak  increa.ses  only  slowly  with  N.  This 
result  implies  that  the  specific  heat  exponent 
a  must  be  very  small  (if  there  is  a  transition) 
or  that  the  correlation  length  exponent  v 
must  be  quite  large.  In  addition,  the  f>eak 
shifts  to  larger  values  of  K  with  increasing 
system  size,  which  is  the  opposite  of  what  is 
usually  observed  for  finite  size  effects  in 
systems  with  periodic  boundary  conditions. 
We  believe  the  most  plausible  explanation  is 
that  the  peak  occurs  when  the  persistence 
length  reaches  the  system  size  and  that 
there  is  no  phase  transition. 

An  additional  confirmation  of  this  inter¬ 
pretation  is  given  by  the  behavior  of  the 
scale-dependent  effective  bending  rigidity 


K,,„.  The  rentirmalized  bending  rigidity  can 
be  expressed  in  terms  of  the  mean  curvature 
susceptibilitv  as  (22) 

\^(0<f/(4)H{^')X  (4) 

where  (A)  is  the  average  membrane  area, 
dent>tcs  the  cuinulant  average,  and  E/(4) 
is  the  mean  curv'ature  at  In  the  ct)ntinuum 
limit,  the  extrinsic  curvature  tensor  is  pro¬ 
portional  to  the  unit  nt)rmal  vector.  The 
di.scretization  (2J)  of  /  d^|V(;(g)i/(^)  we 
use  is  therefore  2;,s.|X„„(CT„//.^.)(r,  -  r,)|, 
where  the  /-sum  runs  over  monomers,  the 
j'(/)-sum  runs  over  the  neighbors  of/,  and  s, 
is  the  sign  of  n,-2,(„(r,  -  r^),  with  n,  the 
surface  normal  at  monomer  /  obtained  by 
averaging  over  the  normals  of  all  triangles 
that  have  /  as  a  corner;  is  the  distance 
between  the  two  neighboring  monomers  /, 
j,  and  <T,y  is  the  length  of  a  bond  in  the  dual 
lattice.  For  the  metric  we  have  chosen, 

=  I/V3  so  that  the  right-hand  side  of  Eq.  4 
becomes 

I  /(()  t  m(k') 

(5) 

For  vesicles,  the  amplitude  of  the  average 
in  Eq.  4  depends  strongly  on  the  si7.e  of  the 
region  over  which  one  averages.  In  particu¬ 
lar,  in  the  continuum  limit,  one  finds  that 
K,.rt  =  6k  when  averaging  over  one-half  of  a 
sphere  when  fluctuations  are  ignored.  Our 
results  for  k  =  k^.h/6,  obtained  by  averaging 
over  N,  =  (N  —  1  )/2  contiguous  sites  on  the 
vesicle  surface,  are  shown  in  Fig.  3B  (24); 
K^.^  is  a  decreasing  function  of  N  for  all  A. 
There  docs  not  appear  to  be  any  qualitative 
difference  in  the  behavior  of  for  large 
and  small  A. 

We  have  not  attempted  to  make  a  quan¬ 
titative  compari-son  with  theory  in  the  large 
A  regime,  because  the  expected  logarithmic 
renormalization  of  k  is  completely  masked 
by  averaging  over  a  finite  portion  of  the 
membrane  surface  (25).  Nevertheless,  the 
essential  qualitative  feature  that  is  a 
decreasing  funaion  of  length  scale  in  both 
the  large  and  small  A  regimes  supports  the 
interpretation  that  there  is  no  phase  transi¬ 
tion  as  a  function  of  A. 

Self-avoiding  fluid  membranes  are  there¬ 
fore  very  likely  always  crumpled  at  large 
length  scales,  independent  of  the  value  of 
the  microscopic  bending  rigidity  k.  Further¬ 
more,  the  crumpled  state  of  the  model  we 
considered  is  a  treelike  ramified  object  char¬ 
acterized  by  ~  N,  in  agreement  with 
results  obtained  for  hypcrcubic  plaquette 
models  {10)  of  random  surfaces.  A  more 
detailed  analysis,  however,  is  required  to 


determine  the  quantitatise  scale  de[X'ndence 
of  the  renormalized  bending  rigiditv. 

Fluid  membranes  with  a  wide  range  of 
bare  bending  rigidities  have  been  studied 
experimentallv  in  considerable  detail  recent- 
ly.  For  phospholipid  bilayers  at  rtxim  tem¬ 
perature,  k/FV  =  10  to  20  (.see  (26-28), 
which  contain  many  additional  references  |, 
so  that  the  persistence  length  is  usually 
much  greater  than  the  size  of  the  membranes 
and  the  reduction  of  the  bending  rigidiU'  by 
shajse  fluctuations  is  negligible.  An  excep¬ 
tionally  small  value  for  however,  is  ob¬ 
tained  for  lecithin  membranes  containing  a 
small  amount  of  bipolar  iipid  (26);  this  may 
al.st)  be  the  case  for  surfactant  bilayers  that 
form  lamellar  and  bicontinuous  phases  (29). 
These  systems  are  therefore  well  suited  for 
experimental  studies  of  flexible  fluid  mem¬ 
branes. 

We  have  shown  that  simulations  per¬ 
formed  using  simple  string-and-bead  mod¬ 
els  for  randomly  triangulated  surfaces  (with 
a  wide  range  of  finite  bare  bending  rigidi¬ 
ties)  can  yield  detailed  information  on  both 
the  large-length-scale  conformation  as  well 
as  the  scale  dependence  of  elastic  constants 
such  as  the  bending  rigidity.  Further  numer¬ 
ical  studies  of  this  model  and  its  generaliza¬ 
tions  should  permit  a  more  detailed  compar¬ 
ison  with  experiment  and  could  provide  new 
insight  into  many  phenomena  involving 
membranes  and  other  self-as.sembling  struc¬ 
tures. 
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